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We address one of the main challenges to Ti02-photocatalysis, namely band gap narrowing, by combining 
nanostructural changes with doping. With this aim we compare Ti02's electronic properties for small OD clus- 
ters, ID nanorods and nanotubes, 2D layers, and 3D surface and bulk phases using different approximations 
within density functional theory and GW calculations. In particular, we propose very small (R < 5 A) but 
surprisingly stable nanotubes with promising properties. The nanotubes are initially formed from Ti02 layers 
with the Pt02 structure, with the smallest (2,2) nanotube relaxing to a rutile nanorod structure. We find that 
quantum confinement effects - as expected - generally lead to a widening of the energy gap. However, substi- 
tutional doping with boron or nitrogen is found to give rise to (meta-)stable structures and the introduction of 
dopant and mid-gap states which effectively reduce the band gap. Boron is seen to always give rise to n-type 
doping while depending on the local bonding geometry, nitrogen may give rise to n-type or p-type doping. For 
undercoordinated Ti02 surface structures found in clusters, nanorods, nanotubes, layers and surfaces nitrogen 
gives rise to acceptor states while for larger clusters and bulk structures donor states are introduced. 



I. INTRODUCTION 

Motivated by the world's ever increasing need for cleaner 
burning fuels and more viable forms of renewable energy, hy- 
drogen production via photocatalysis has been intensely re- 
searched as a possible candidate for addressing these issues. 
Since the first experimental formation of hydrogen by pho- 
tocatalysis in the early 1980s [1], Ti02 has been the catalyst 
of choice. Reasons for this include the position of Ti02's con- 
duction band above the energy of hydrogen formation, the rel- 
atively long lifetime of excited electrons which allows them to 
reach the surface from the bulk, Ti02's high corrosion resis- 
tance compared to other metal oxides, and its relatively low 
cost [2-4]. 

However, the large band gap of bulk Ti02 (~ 3 eV) means 
that only high energy UV light may excite its electrons. This 
effectively blocks most of the photons which pierce the atmo- 
sphere, typically in the visible range, from participating in any 
bulk Ti02 based photocatalytic reaction. On the other hand, 
the difference in energy between excited electrons and holes, 
i.e. the band gap, must be large enough (> 1.23 eV) to dis- 
sociate water into hydrogen and oxygen. For these reasons it 
is of great interest to adjust the band gap £gap of Ti02 into 
the range 1.23 < Sgap ^ 2.5 eV, while maintaining the useful 
properties mentioned above [5]. 

With this aim, much research has been done on the influ- 
ence of Ti02 nanostructure [6-10] and dopants [5, 11-17] on 
photocatalytic activity. For low dimensional nano structured 
materials, electrons and holes have to travel shorter distances 
to reach the surface, allowing for a shorter quasi-particle life- 
time. However, due to quantum confinement effects, lower di- 
mensional Ti02 nanostructures tend to have larger band gaps 
[18]. On the other hand, although doping may introduce mid- 
gap states, recent experimental studies have shown that boron 
and nitrogen doping of bulk Ti02 yields band gaps smaller 
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than the threshold for water splitting [11, 12]. This suggests 
that low dimensional structures with band gaps larger than 
about 3.0 eV may be a better starting point for doping. 

The experimental synthesis and characterization of nanos- 
tructured materials is in general a costly and difficult task. 
On the other hand, modern electronic structure modelling has 
reached a level where large-scale calculations can provide re- 
alistic descriptions of structure and electronic properties. 

Based on our investigation, we suggest as a promising can- 
didate small {R < 5 k) T1O2 nanotubes, with a hexagonal 
ABC Pt02 structure (HexABC), which we find to be surpris- 
ingly stable, even in the boron and nitrogen doped forms. This 
stability may be attributed to their structural similarity to bulk 
rutile Ti02, with the smallest nanotube having the same struc- 
ture as a rutile nanorod. 

A further difficulty for any photocatalytic system is control- 
ling how electrons and holes travel through the system [19]. 
For this reason, methods for reliably producing both n-type 
and p-type Ti02 semiconducting materials are highly desir- 
able. So far, doped Ti02 tends to yield only n-type semicon- 
ductors. We propose that p-type Ti02 semiconducting mate- 
rials may be obtained by nitrogen doping surface sites of low 
dimensional materials. 

In this study we report the results of density functional 
theory (DFT) and GW calculations of the energetic stability 
and electronic structure of recently suggested (Ti02)n clus- 
ters (n < 9) [6] and novel (n,n) Ti02 nanotubes (n < 4) in 
the undoped, boron doped, and nitrogen doped forms. The 
formation energy Eform, density of states (DOS), and energy 
gap Sgap, for these systems are compared with that for 2D 
HexABC and anatase layers, and 3D Ti02 anatase surface, 
anatase bulk, and rutile bulk phases. Furthermore, we analyze 
how boron and nitrogen doping influences the DOS for these 
systems, and how their nanostructure may determine whether 
the resulting semiconductor is n-type or p-type. 

In Sec. II we describe the DFT and GW methodologies used 
to obtain the energies and electronic structure of the studied 
systems. We compare the energetic stability of the systems 
in Sec. III. In Sec. IV we discuss the electronic structures of 
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the systems, showing the DOS and energy gaps, followed by 
a concluding section. 



II. METHODOLOGY 

All DFT calculations have employed the RPBE exchange 
correlation (xc)-functional [20]. The plane-wave code DA- 
CAPO [21, 22] was used for structural minimization, the 
real-space code OCTOPUS [23] for charged calculations, and 
YAMBO[24] with ABINIT[25] or PWSCF[26] for GW calcu- 
lations. A plane- wave cutoff of 340 eV was used, with a 
Monkhorst-Pack k-point sampling of 1 x 1 x 12 for Ti02 nan- 
otubes, where the nanotube axis is parallel to the z-axis, 
12x12x1 for Ti02 layers and surfaces, where the normal 
direction is parallel to the z-axis, and 12x12x12 for Ti02 
bulk phases. All structures have been relaxed until a maxi- 
mum force below 0.04 eV/A was obtained. The occupation 
of the one electron states was calculated at a temperature of 
ksT ^0.1 eV for the periodic systems and ksT ^0.01 eV 
for the clusters, with all energies extrapolated to T =0 K. 
Spin unpolarized calculations have been performed for the un- 
doped Ti02 systems, while spin polarized calculations have 
been performed for all doped Ti02 systems, since the unit 
cells for the doped systems contain an odd number of elec- 
trons. 

Doping of (Ti02)n clusters has been modeled by substitut- 
ing a single boron or nitrogen atom in each geometrically in- 
equivalent oxygen site of the most stable isomer to determine 
the most stable doping site. Only clusters of sufficient size 
to obtain experimentally realizable doping fractions < 10% 
(5 < n < 9) have been considered [11-15]. 

To model Ti02 nanotube doping, we have repeated the min- 
imal unit cell four times along the tube axis, and substituted a 
single boron or nitrogen atom in each geometrically inequiv- 
alent oxygen site to obtain the most stable doped structure. In 
this way, dopant-dopant interactions are minimized by main- 
taining a dopant separation of approximately 12 A. This corre- 
sponds to experimentally realizable doping fractions of 3.1%, 
2.1%, and 1.6% for Ti02 (2,2) nanorods, (3,3) nanotubes and 
(4,4) nanotubes, respectively. 

Doping of Ti02 layers, surfaces, and bulks has been simi- 
larly modeled by repeating the minimal unit cell twice in each 
periodic direction, and substituting a single boron or nitrogen 
atom in each geometrically inequivalent oxygen site to obtain 
the most stable structure. Experimentally realizable doping 
fractions of 5.6%, 3.1%, and 3.1% were thus obtained for the 
layers, surfaces, and bulks respectively. 
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FIG. 1: Formation energy Eform in eV per Ti02 functional unit 
versus Ti02 structure for OD (Ti02)n clusters (n < 9), ID Ti02 
(2,2) nanorods, (3,3) nanotubes, and (4,4) nanotubes, 2D HexABC 
and anatase layers, and 3D anatase surface, anatase bulk, and rutile 
bulk phases. DFT calculations using RPBE for the undoped (O, red), 
nitrogen doped (□, blue), and boron doped (A, magenta) systems are 
shown. Note the sp^ character of oxygen in the more stable rutile 
and (2,2) nanorod structure, and the sp^ character of oxygen in the 
less stable layer and larger nanotube structures. 



III. ENERGETIC STABILITY 

We define the formation energy Eform for a given structure 
consisting of n Ti02 functional units as 

Eform = -^[Tin02n-p-gBpNg] - ^^[TiOa] 

-- {pE[B] + qE[N] - (p + ^)^[0]) , (1) 



where E[Tin02n-p-q^p^q] is the total energy for the system, 
£^[Ti02] is the energy of an isolated Ti02 molecule, and E[0], 
E[B], and E[N] are the respective energies with reference to 
gas phase species for O, B, and N. These have been obtained 
using the experimental doping reactions 

^[O] = ^[H20] - E[ll2] - Ai7[H20], (2) 
^[B] = i(^[B2H6]-3^[H2]-Ai7[B2H3]), (3) 
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where the formation reaction enthalpies are Ai^[H20] ^ 
-2.506 eV and Ai^[B2H6] ^ 0.377 eV, as taken from 
Ref. [27]. This avoids difficulties associated with modeling 
isolated atoms and the triplet state of molecular oxygen. 

Fig. 1 shows DFT calculated formation energies Eform 
for undoped, boron doped, and nitrogen doped forms of OD 
(Ti02)n clusters (n < 9), ID (n,n) nanotubes (n < 4), 2D 
Hex ABC and anatase layers, and 3D anatase surface, anatase 
bulk, and rutile bulk phases. 

From the formation energies, we see that (Ti02)n clusters 
become generally more stable with increasing size, as ex- 
pected. On the other hand, the formation energies of the (n,n) 
Ti02 nanotubes (n < 4) suggest that they are all surprisingly 
stable, being more stable than the 2D structures and within ap- 
proximately 0.2 eV of the least stable bulk phases. Even more 
surprising, the Ti02 nanotubes seem to become more stable 
with decreasing size. This is clearly seen in Fig. 2, which 
shows the radial dependence of the formation energies for the 
Ti02 (2,2) nanorod and Ti02 (n, n) nanotubes (3 < n < 8). 

In a simple classical picture, the strain energy to be over- 
come when bending a layer into a tube should be inversely 
proportional to the square of the tube radius (~ R~'^). Such 
models have been shown to accurately describe Ti02 anatase 
and lepidocrocite layer nanotubes [28]. However, the Hex- 
ABC Ti02 nanotube system is complicated by the fact that 
Ti02 is not found naturally in a Hex ABC layered structure. 
As a result, the strain energy introduced by bending the layer 
may be overcome by changes in bonding which better resem- 
ble the bulk Ti02 structures. This is seen in Fig. 1, where we 
find oxygen atoms have sp^ character in both the infinite layer 
structure and the (4,4) Ti02 nanotube, but sp^ character in the 
more stable bulk rutile and (2,2) Ti02 nanorod. 

As shown in Fig. 2, the radial dependence of the nanotube 
formation energies may be well approximated by 
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Here E layer ^4.63 eV/Ti02 is the formation energy of the 
nanotube layer calculated using the optimized unit cell pa- 
rameters for the nanotubes, EmtUe ^ 4.89 eV/Ti02 is the 
formation energy of the bulk rutile phase, R(2,2) the radius 
of a Ti02 (2,2) nanorod, and R is the nanotube radius. The 
first term in (5) ensures the nanotube formation energy tends 
to that for the infinite layer at large nanotube radii. The sec- 
ond term models the increase in stability as the nanotube is 
'bent' into a more rutile-like structure at smaller radii. This 
reaches a maximum for a (2,2) nanorod (R = R{2,2))^ where 
the structure becomes that of the rutile phase. 

It should be noted that DFT formation energies have been 
shown to not accurately describe the relative stability of 
anatase and rutile phases, independent of xc-functional used 
[29] . This should be borne in mind when comparing forma- 
tion energies. 

Comparison of the doped formation energies with those of 
the undoped Ti02 structures shows that although less stable, 
doped forms of both Ti02 clusters and nanotubes should be 
experimentally realizable. This is further justified by recent 
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FIG. 2: Formation energy Eform in eV per Ti02 functional unit 
versus radius Rin A for Ti02 (2,2) nanorods and (n, n) nanotubes 
(3 < n < 8) (O, red). For comparison, the Ti02 nanotube formation 
energies are also approximated by Eform ~ Eiayer + [ErutUe - 

Eiayer)R{2,2)/ R, ( , red), whcrc R{2,2) ~2.54 A is the radius of 

a Ti02 (2,2) nanorod. 



nitrogen doping experiments of large Ti02 nanotubes which 
yielded doping fractions between 1% and 10% [12, 14, 15]. 

The formation energies of the doped systems also appear 
to correlate well with both those of the undoped structures, 
and the relative doping fraction. Further, structures are typ- 
ically less stable when doping induces oxygen dislocations, 
which occur for boron and nitrogen doped (Ti02)7 clusters, 
(Ti02)8 clusters, and bulk phases. On the other hand, systems 
are more stable when the Ti02 structure is unchanged by the 
introduction of a dopant, as is the case for nitrogen doping of 
(Ti02)n clusters with surface sites in = 5,6,9), nanotubes, 
layers, and surfaces. Such a distinction will also prove useful 
when we analyze the electronic structure of the doped systems 
in Sec. IV. 



IV. ELECTRONIC STRUCTURE 
A. DOS of OD Clusters 

We shall begin our analysis of the electronic structure of 
boron and nitrogen doped Ti02 nanostructures by first con- 
sidering (Ti02)n clusters (5< n <9). As well as being of 
interest in its own right, the study of such systems provides 
insight into how local changes in nanostructure may change 
whether a dopant behaves as a donor or acceptor site. Such in- 
formation will prove useful when trying to form p-type Ti02 
semiconducting materials. 

The DFT calculated DOS and structures for the most sta- 
ble boron and nitrogen doped (Ti02)n clusters are shown in 
Figs. 3 and 4, for 5 < n < 9. In general, boron prefers to 
replace the most highly coordinated oxygen in Ti02 clusters. 
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FIG. 3: Total density of states in eV ' per TiOa functional unit versus 
energy in eV for undoped (thin black solid line), boron doped (ma- 
genta dashed line) and nitrogen doped (blue thick solid line) (a-e) 
(Ti02)n clusters for 5 < n < 9, with doping fractions of i/2n. Ener- 
gies are measured from the top of the valence bands of anatase Ti02 
svB, and the DOSs are shifted to align the lowest eigenstate with that 
for anatase Ti02. Occupancy is denoted by curve filling for states in 
the band gap. 
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FIG. 4: Schematics of undoped (left), boron doped (middle), and ni- 
trogen doped (right) (Ti02)n clusters for 5 < n < 9, with doping 
fractions of ^/2n. Note the circled nitrogen-oxygen bonds in the ni- 
trogen doped structures for n =7 and 8, where nitrogen acts as a 
donor [6]. 



typically in a central location, and forms boron-oxygen bonds 
via oxygen dislocations, as shown in Fig. 4. Due to its elec- 
tropositive character, boron acts as a donor, introducing three 
occupied mid-gap states near the lowest unoccupied molecu- 
lar orbital (LUMO), as seen in Fig. 3. 

However, for (Ti02)n clusters the influence of nitrogen 
dopants is not as straightforward. As with rutile Ti02 
surfaces [16, 17], nitrogen prefers to occupy sites which are 
3 -fold coordinated to titanium. For clusters consisting only 
of surface sites (n = 5, 6, 9), oxygen is tightly constrained, 
and no significant changes in the cluster's structure occur via 
oxygen dislocations, as shown in Fig. 4. As with rutile Ti02 
surfaces [16, 17], nitrogen acts as an acceptor, introducing an 
unoccupied mid-gap state near the highest occupied molecular 
orbital (HOMO). 



For clusters where the nitrogen dopant's oxygen neighbors 
occupy interior/bulk sites (n = 7, 8), oxygen is more mo- 
bile, and may form new nitrogen-oxygen bonds, as shown in 
Fig. 4. Since oxygen is more electronegative than nitrogen, 
nitrogen may transfer charge to the bonding oxygen through 
these bonds. In this case, nitrogen instead acts as a donor, 
introducing occupied states near the LUMO. 

In general, we find when dopants form bonds with oxy- 
gen atoms via oxygen dislocations, they act as donor sites. 
On the other hand, when dopants such as nitrogen leave the 
Ti02 structure unchanged, they act as acceptor sites, yield- 
ing a somewhat more stable structure. We shall find such a 
description useful in understanding why nitrogen dopants be- 
have differently in low dimensional systems and bulk systems. 
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FIG. 5: RPBE calculation of the total density of states in eV" per 
Ti02 functional unit versus energy in eV for undoped (thin black 
solid line), boron doped (magenta dashed line) and nitrogen doped 
(blue thick solid line) Ti02 (a) (2,2) nanorods, (b) (3,3) nanotubes 
and (c) (4,4) nanotubes, with 3.1%, 2.1%, and 1.6% doping respec- 
tively. Energies are measured from the top of the valence bands of 
anatase Ti02 svb, and the DOSs are shifted to align the lowest eigen- 
state with that for anatase Ti02. Occupancy is denoted by curve fill- 
ing for states in the band gap. 




(3,3) Nanotube 




(4,4) Nanotube 




Q B-doped 



N-doped 



FIG. 6: Schematics of boron doped (left) and nitrogen doped (right) 
Ti02 (2,2) nanorods, (3,3) nanotubes, and (4,4) nanotubes, with 
3.1%, 2.1%, and 1.6% doping respectively. The highest occupied 
states for boron and nitrogen doped Ti02 ID structures are depicted 
by isosurfaces of ±0.05 e/A^. 



B. DOS of ID Nanorods and Nanotubes 

We shall now discuss how boron and nitrogen doping in- 
fluences the electronic structure of stable ID nanorods and 
nanotubes. Although the the most stable doping sites are the 
same as those obtained for doped Ti02 rutile surface and bulk 
[5, 16, 17], the influence of dopants on the DOS is rather dif- 
ferent. 

Figure 5 shows the DFT calculated DOS and Fig. 6 the 
structures of the most stable boron doped and nitrogen doped 
Ti02 (2,2) nanorods, (3,3) nanotubes, and (4,4) nanotubes. 
The highest occupied state is also shown as isosurfaces of 
±0.05e/A^ in the side views of the doped structures. 

As with Ti02 clusters, the influence of boron dopants on 
Ti02 nanorods and nanotubes may be understood in terms 
of boron's weak electronegativity, especially when compared 
with the strongly electronegative oxygen. We find that boron 
prefers to occupy oxygen sites which are 2-fold coordinated 
to neighboring titanium atoms. However, as with the OD clus- 
ters, boron's relatively electropositive character induces sig- 
nificant structural changes in the ID structures, creating a 
stronger third bond to a neighboring three-fold coordinated 
oxygen via an oxygen dislocation, as shown in Fig. 5. This 
yields three occupied mid-gap states localized on the boron 
dopant, which overlap both the valence band O 2p-j^ and con- 
duction band Ti 2dxy states, as shown in Fig. 5. Boron dopants 




FIG. 7: GW calculation of the total density of states in eV"^ per Ti02 
functional unit versus energy in eV for undoped (thin black solid 
line), boron doped (magenta dashed line) and nitrogen doped (blue 
thick solid line) Ti02 (2,2) nanorods with 3.1% doping. Energies 
are measured from the top of the valence bands of the undoped Ti02 
nanorod svb, and the DOSs are shifted to align the lowest eigenstate 
with that for the undoped Ti02 nanorod. Occupancy is denoted by 
curve filling for states in the band gap. 
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FIG. 8: Total density of states in eV"^ per Ti02 functional unit ver- 
sus energy in eV for undoped (thin black solid line), boron doped 
(magenta dashed line) and nitrogen doped (blue thick solid line) (a) 
Ti02 HexABC layer, (b) Ti02 anatase surface, and (c) bulk anatase, 
with 5.6%, 3.1%, and 3.1% doping respectively. Energies are mea- 
sured from the top of the valence bands of anatase Ti02 svb, and the 
DOSs are shifted to align the lowest eigenstate with that for anatase 
Ti02. Occupancy is denoted by curve filling for states in the band 
gap. 



thus yield donor states near the conduction band, which may 
be photocatalytically active in the visible region. However, 
the quantum confinement inherent in these ID structures may 
stretch these gaps, as found for the GW calculated DOS shown 
in Fig. 7. 

On the other hand, we find nitrogen dopants prefer to oc- 
cupy oxygen sites which are 3 -fold coordinated to Ti, as was 
previously found for the rutile Ti02 surface[16, 17]. This 
yields one occupied state at the top of the valence band and 
one unoccupied mid-gap state in the same spin channel. Both 
states are localized on the nitrogen dopant but overlap the va- 
lence band O 2^^^ states, as shown in Fig. 6. Nitrogen dopants 
thus act as acceptors, providing localized states well above 
the valence band, as is also found for the GW calculated DOS 
shown in Fig. 7. 

Although we find nitrogen dopants act as acceptors in Ti02 
ID structures, such large gaps between the valence band and 
the unoccupied mid-gap states would not yield p-type semi- 
conductors. This may be attributed to the substantial quantum 
confinement in these ID structures. However, for 2D and 3D 
systems, we find it is possible to produce both p-type and n- 
type classical semiconductors. 
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FIG. 9: Schematics of the boron doped (left) and nitrogen doped 
(right) HexABC layer, anatase surface, and anatase bulk phases, with 
5.6%, 3.1%, and 3.1% doping respectively. The highest occupied 
states for boron and nitrogen doped Ti02 2D and 3D structures are 
depicted by isosurfaces of =b0.05e/A^. Note the circled nitrogen- 
oxygen bonds in the nitrogen doped structures for the bulk phase, 
where nitrogen acts as a donor [6]. 



C. DOS of 2D Layers, 3D Surfaces and Bulks 

To determine the reliability of our DFT and GW electronic 
structure calculations, we now consider the influence of boron 
and nitrogen doping on the DOS for 2D and 3D Ti02 struc- 
tures, which may be compared with recent experiments. 

For the HexABC layer shown in Figs. 8(a) and 9, boron 
and nitrogen dopants act as donors and acceptors respectively, 
as was the case for the Ti02 nanotubes. However, for these 
2D structures the doping states are sufficiently close in en- 
ergy to the conduction band or valence band to allow charge 
transfer for classical n-type or p-type doping by boron or ni- 
trogen dopants respectively. A similar behavior is found for 
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the anatase surface, shown in Fig. 8(b) and Fig. 9. This sug- 
gests boron and nitrogen doped 2D Ti02 layers or surfaces 
may potentially be used for Ti02 based electronics. 

On the other hand, for bulk anatase boron and nitrogen 
dopants both yield n-type semiconductors, as is seen exper- 
imentally [13]. We attribute this difference to the nitrogen- 
oxygen bond shown in Fig. 8(c), which may arise due to oxy- 
gen's greater structural flexibility in the bulk. As a result, 
boron and nitrogen doping both yield mid-gap occupied Ti 
2dxy conduction band states with little weight on the doping 
sites. Here, the dopants donate their valence electrons fully to 
the conduction band, as seen in Fig. 8(c). 

Nitrogen thus acts as a donor for bulk-like clusters and 
anatase, in agreement with recent findings for nitrogen 
doping of both bulk anatase and thick bulk-like anatase 
nanopores [13]. This suggests that n-type or p-type semi- 
conductors may be produced by nitrogen doping depending 
on whether nearby oxygen atoms occupy lower-coordinated 
surface or higher-coordinated bulk sites. 

Further, we also find that the highest occupied states for 
boron and nitrogen doped systems become increasingly local- 
ized as the structure is dimensionally constrained from bulk 
surface layer nanotube nanorod. This may have im- 
portant consequences for the photocatalytic activity of these 
more localized states in lower dimensional structures. On the 
other hand, this may be partially alleviated by having these 
states located on the structure's surface. In this way, electron- 
hole pairs need not travel significantly, as is the case for bulk 
materials. 



D. Energy Gaps 

To analyze the qualitative trends in photocatalytic activity 
of boron and nitrogen doped Ti02 nanostructures, we shall 
now compare the energy gaps obtained using different approx- 
imations within DFT, GW calculations, and experiment. 

The DFT calculated energy gaps Sgap, for (Ti02)n clusters 
(n < 9) and band gaps for Ti02 (2,2) nanorods, (3,3) nan- 
otubes, and (4,4) nanotubes, in the undoped, boron doped, and 
nitrogen doped forms, are shown in Fig. 10. For the clusters, 
the gap is estimated by both the difference in energy between 
the HOMO and LUMO, and the difference between the ion- 
ization potential and electron affinity energies Ip — Ea. The 
band gap for Ti02 nanorods, nanotubes, layers, surfaces and 
bulk phases is approximated by the indirect band gap between 
highest occupied and lowest unoccupied states. Although not 
as relevant for photoabsorption as the direct gap, the large 
size required for the doped super cells, shown in Fig. 6 and 
Fig. 9, leaves indirect and direct gaps indistinguishable. For 
the undoped nanorods and nanotubes and bulk rutile we find 
the band gap is direct, while for the undoped Hex ABC layer 
and anatase structures we find the band gap is indirect. For 
the same structures we also provide B3LYP band gaps, GW 
calculations, and experimental results for comparison. 

As previous studies have shown, standard DFT tends to 
underestimate band gaps for bulk Ti02 by approximately 1 
eV, due in part to self-interaction errors [31, 32]. This may 
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FIG. 10: Energy gap Sgap in eV versus Ti02 structure for OD (Ti02)n 
clusters (n < 9), ID Ti02 (2,2) nanorods, (3,3) nanotubes, (4,4) 
nanotubes, 2D HexABC and anatase layers, and 3D anatase surface, 
anatase bulk, and rutile bulk phases. DFT calculations using RPBE 
of the highest occupied and lowest unoccupied state gaps / Ip — Ea 
for the (a) undoped (0/#, red), (b) boron doped (A/A, magenta), and 
(c) nitrogen doped (□/■, blue) systems are compared with undoped 
B3LYP (T), GW (►), and experimental (u) results [6, 12, 18, 30]. 
Small open symbols denote transitions between highest fully occu- 
pied states and the conduction band. 



be partially addressed by the use of hybrid functionals such 
as B3LYP, which generally seem to improve band gaps for 
bulk systems [33-35]. However, B3LYP calculations for Ti02 
clusters largely underestimate the gap relative to the more re- 
liable Ip — Ea, while for the bulks the gap is overestimated by 
approximately 0.4 eV relative to experiment. Also, B3LYP 
and RPBE calculations provide the same qualitative descrip- 
tion of the trends in the energy gaps for Ti02. 

GW is probably the most successful and generally appli- 
cable method to calculate quasi-particle gaps. For clusters it 
agrees well with Ip — Ea, but for bulk systems GW seems to 
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overestimate the experimental gap like B3LYP. This overesti- 
mation by GW may be attributed to excitonic effects, which 
are not included in GW, or the role played by oxygen defects 
[36]. For these reasons, GW and DFT calculations are best 
used as upper and lower bounds for optical band gaps. 

Fig. 10 shows that for both 3D and 2D systems, RPBE gaps 
underestimate the experimental results by approximately 1 eV. 
For ID and OD systems, we find a much larger difference of 
about 4 eV and 5 eV respectively, between the RPBE gaps and 
the Ip — Ea and GW results. We may attribute this increas- 
ing disparity to the greater quantum confinement and charge 
localization in the ID and OD systems, which yield greater 
self-interaction effects. We find B3LYP gaps also underesti- 
mate this effect, simply increasing the RPBE energy gaps for 
both OD and 3D systems up by about 1 .4 eV. 

On the other hand, we find the RPBE gaps reproduce qual- 
itatively the structural dependence of the Ip — Ea, GW, and 
experimental results for a given dimensionality, up to a con- 
stant shift. This is true even for 3D bulk systems, where stan- 
dard DFT does not predict rutile to be the most stable [29], as 
shown in Fig. 1. 

Whether calculated using RPBE, Ip — Ea or GW, the en- 
ergy gaps for both boron and nitrogen doped Ti02 nanostruc- 
tures are generally narrowed, as shown in Fig. 10(b) and (c). 
However, for nitrogen doped (Ti02)n clusters where nitrogen 
acts as an acceptor (n = 5, 6, 9), the energy gap is actually 
increased when spin is conserved, compared to the undoped 
clusters in RPBE. This effect is not properly described by the 

^ + 1 transitions of Ip — Ea, for which spin is not con- 
served for these nitrogen doped clusters. On the other hand, 
when nitrogen acts as a donor (n = 7, 8) the smallest gap 
between energy levels does conserve spin. 

The boron doped Ti02 nanorods and nanotubes have per- 
haps the most promising energy gap results of the Ti02 struc- 
tures considered herein, as seen in Fig. 5(b). Boron dopants 
introduce in the nanorods localized occupied states near the 
conduction band edge in both RPBE (c/ Figs. 5 and 6) and 
GW (c/ Fig. 7) calculations. On the other hand, nitrogen dop- 
ing of nanorods and nanotubes introduces well defined mid- 
gap states, as shown in Fig. 5(c). However, to perform water 
dissociation, the energy of the excited electron must be above 
that for hydrogen evolution, with respect to the vacuum level. 
This is not the case for such a mid-gap state. This opens the 
possibility of a second excitation from the mid-gap state to 
the conduction band. However, the cross section for such an 



excitation may be rather low. 

For boron doping of 2D and 3D structures, the highest oc- 
cupied state donates its electron almost entirely to the con- 
duction band, yielding an n-type semiconductor. Thus at very 
low temperatures, the RPBE band gap is very small. The same 
is true for n-type nitrogen doped bulk anatase. For these rea- 
sons, we also provide the energy gap between the highest fully 
occupied state and the conduction band, which may be more 
relevant for photoabsorption. We find these RPBE gaps are 
still generally smaller than those for their undoped Ti02 coun- 
terparts. 

In summary, for both boron and nitrogen doped clusters we 
find RPBE gaps differ from Ip — Ea by about 3 eV, while for 
nitrogen doped anatase the RPBE gap differs from experiment 
by about 0.6 eV. Given the common shift of 1 eV for undoped 
2D and 3D structures, this suggests that both boron and ni- 
trogen doped 2D Ti02 structures are promising candidates for 
photocatalysis. Further, the boron and nitrogen doped ID nan- 
otube results also warrant further experimental investigation. 



V. CONCLUSION 

In conclusion, we have demonstrated how the electronic 
properties of Ti02 may be "tailored" using nanostructural 
changes in combination with boron and nitrogen doping. 
While boron doping tends to produce smaller band gap n- 
type semiconductors, nitrogen doping produces p-type or n- 
type semiconductors depending on whether or not nearby oxy- 
gen atoms occupy surface sites. This suggests that a p-type 
Ti02 semiconductor may be produced using nitrogen doping 
in conjunction with surface confinement at the nanoscale. 
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